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The Color Glass Condensate 3 
1 High energy QCD and the Color Glass Condensate 



Only 5% of the mass of the universe is "bright matter" , but 99% of this visible 
matter is described by QCD, the theory of the strong interactions. QCD has 
been described as a nearly perfect theory, with the only parameters being the 
masses of the current quarks (1). The vast bulk of visible matter is therefore 
"emergent" phenomena arising from the rich dynamics of the QCD vacuum and 
the interactions of the fundamental quark and gluon constituents of QCD. Al- 
though we have enough information from a wide range of experimental data and 
from numerical lattice computations that QCD is the right theory of the strong 
interactions, outstanding questions remain about how a variety of striking phe- 
nomena arise in the theory. This is because the complexity of the theory also 
makes it very difficult to solve. 

Our focus in this review is on high energy scattering in QCD. A traditional 
approach to these phenomena is to divide them into "hard" or "soft" scatter- 
ing, corresponding respectively to large or small momentum exchanges in the 
scattering. In the former case, because of the "asymptotic freedom" of QCD, 
phenomena such as jets can be computed in a perturbative framework. In the 
latter case, because the momentum transfer is small, the "infrared slavery" of 
QCD suggests that the coupling is large; the scattering therefore is intrinsically 
non-perturbative and therefore not amenable to first principles analysis. This is 
problematic because "soft" dynamics comprises the bulk of QCD cross-sections. 
In contrast, the perturbatively calculable hard cross-sections are rare processes 1 . 

We shall argue here that the traditional separation of hard versus soft QCD 
dynamics is oversimplified because novel semi-hard scales generated dynamically 
at high energies allow one to understand highly non-perturbative phenomena in 
QCD using weak coupling methods. To clarify what we mean, consider inclusive 
cross-sections in deeply inelastic scattering (DIS) experiments of electrons off 
hadrons. 

Inclusive cross-sections in DIS (see figure 1) can be expressed in terms of the 
Lorentz invariants i) x which corresponds, at lowest order in perturbation theory, 
to the longitudinal momentum fraction carried by a parton in the hadron, ii) the 
virtual photon four-momentum squared q 2 = —Q 2 < exchanged between the 
electron and the hadron, hi) the inelasticity y, the ratio of the photon energy to 

lattice QCD is not of much help here because it is best suited to compute static properties 

of the theory. 




Figure 1: Kinematics of Deep Inelastic Scattering. 

the electron energy in the hadron rest frame, and iv) the center of mass energy 
squared s. These satisfy the relation x = Q 2 /sy. Instead of the dichotomy of 
"hard" versus "soft", it is useful, for fixed y, to consider two asymptotic limits in 
DIS that better illustrate the QCD dynamics of high energy hadron wavefunc- 
tions. The first, called the Bjorken limit, corresponds to fixed x with Q 2 , s — > oo. 
The second is the Regge-Gribov limit of fixed Q 2 , x — > and s — > oo. 

In the Bjorken limit of QCD, one obtains the parton model wherein the hadron 
is viewed, in the infinite momentum frame (IMF), as a dilute collection of va- 
lence quarks and "wee" (in the terminology coined by Feynman) partons-small 
x gluons and sea quark pairs. Albeit the number of wee partons is large at high 
energies, the hadron is dilute (as illustrated in figure 2) because the phase space 
density is very small for Q 2 — > oo. In this limit, to leading order in the coupling 
constant, the interaction of a probe with the hadron, on the characteristic time 
scale 1/Q, can be expressed as a hard interaction with an individual parton. In 
this "impulse approximation" , the interaction of the struck parton in the hadron 
with co-moving partons is suppressed due to time dilation. The separation of 
hard and soft scattering alluded to previously is valid here. Powerful tools such 
as the Operator Product Expansion (OPE) and factorization theorems extend, to 
high orders in the coupling constant expansion, the hard-soft separation between 
process-dependent physics at the scale 1/Q and universal, long distance, non- 
perturbative Parton Distribution Functions (PDFs). The evolution of the sepa- 
ration of hard and soft scales is given by renormalization group (RG) equations 
called the DGLAP (Dokshitzer-Gribov-Lipatov-Altarelli-Parisi) equations (2). 

The study of strong interactions in the Regge-Gribov limit predates QCD, and 
underlies concepts such as soft Pomeron and Reggeon exchanges, whose dynam- 
ics, described by "Reggeon Field Theory" (RFT), was believed to encompass 
much of the phenomena of multi-particle production. Albeit phenomenologically 




Figure 2: Left: the x-evolution of the gluon, sea quark, and valence quark distri- 
butions for Q 2 = 10 GeV 2 measured at HERA (3). Right: the "phase-diagram" 
for QCD evolution; each colored dot represents a parton with transverse area 
5S± ~ 1/Q 2 and longitudinal momentum k + = xP + . 

very suggestive (4), the dynamics of RFT is intrinsically non-perturbative and 
therefore not easily amenable to systematic computations. With the advent of 
high energy colliders, hadron structure in the Regge-Gribov limit can be explored 
with Q 2 > 1 GeV 2 . As strikingly demonstrated by the HERA DIS data shown in 
figure 2, the gluon distribution xG(x, Q 2 ) in a proton rises very fast with decreas- 
ing x at large, fixed Q 2 - roughly, as a power x~ x with A ~ 0.3. In the IMF frame 
of the parton model, xG(x,Q 2 ) is the number of gluons with a transverse area 
SS ± > l/Q 2 and a fraction k+/P+ ~ x of the proton longitudinal momentum . 
In the Regge-Gribov limit, the rapid rise of the gluon distribution at small x is 
given by the BFKL (Balitsky-Fadin-Kuraev-Lipatov) equation (5), which we will 
discuss at some length later. 

The stability of the theory formulated in the IMF requires that gluons have 
a maximal occupation number of order l/ct s . This bound is saturated for gluon 
modes with transverse momenta k± < Q s , where Q s (x) is a semi-hard scale, 
the "saturation scale", that grows as x decreases. In this novel "saturation" 
regime of QCD (6), illustrated in figure 2 (right panel), the proton becomes a 
dense many body system of gluons. In addition to the strong x dependence, 
the saturation scale Q s has an A dependence because of the Lorentz contraction 
of the nuclear parton density in the probe rest frame. The dynamics of gluons 
in the saturation regime is non-perturbative as is typical of strongly correlated 
systems. However, in a fundamental departure from RFT, this dynamics can be 
The light cone co-ordinates are defined as fc* = (fc° ± fe 3 )/\/2. 
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computed using weak coupling methods as a consequence of the large saturation 
scale dynamically generated by gluon interactions. Thus instead of the "hard" 
plus "soft" paradigm of the Bjorken limit, one has a powerful new paradigm in 
the Regge-Gribov limit to compute the bulk of previously considered intractable 
scattering dynamics in hadrons and nuclei. 

The Color Glass Condensate (CGC) is the description of the properties of 
saturated gluons in the IMF in the Regge-Gribov limit. The effective degrees 
of freedom in this framework are color sources p at large x and gauge fields 
at small x. At high energies, because of time dilation, the former are frozen 
color configurations on the natural time scales of the strong interactions and 
are distributed randomly from event to event. The latter are dynamical fields 
coupled to the static color sources. It is the stochastic nature of the sources, 
combined with the separation of time scales, that justify the "glass" appellation. 
The "condensate" designation comes from the fact that in the IMF, saturated 
gluons have large occupation numbers 0(l/a s ), with typical momenta peaked 
about a characteristic value k± ~ Q s . The dynamical features of the CGC are 
captured by the JIMWLK 3 renormalization group (RG) equation that describes 
how the statistical distribution W[p] of the fast sources at a given x scale evolves 
with decreasing x. The JIMWLK RG is Wilsonian in nature because weakly 
coupled fields integrated out at a given step in the evolution are interpreted as 
"induced color charges" that modify the statistical weight distribution. 

The CGC framework is quite powerful because, given an initial non-pertur- 
bative distribution of sources at an initial scale xq, it allows one to compute 
systematically n-point gluon correlation functions and their evolution with x 
order by order in perturbation theory. In analogy to parton distribution functions 
in the Bjorken limit, the distribution W[p] captures the properties of saturated 
gluons. Unlike PDFs however, which in the IMF are parton densities, W[p] 
carries much more information. Further, in contrast to the "twist" expansion 
of the Bjorken limit which becomes extremely cumbersome beyond the leading 
order in l/Q, the CGC framework includes all twist correlations. Like the PDFs, 
the W's are universal: the same distributions appear in computations of inclusive 
quantities in both lepton-nucleus and hadron-nucleus collisions. 

The CGC is interesting in its own right in what it reveals about the collective 
dynamics of QCD at high parton densities. The CGC RG equations indicate 
that -at fixed impact parameter- a proton and a heavy nucleus become indis- 

3 JIMWLK = Jalilian-Marian, Iancu, McLerran, Weigert, Leonidov, Kovner. 
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tinguishable at high energy; the physics of saturated gluons is universal and 
independent of the details of the fragmentation region. This universal dynam- 
ics has a correspondence with reaction-diffusion processes in statistical physics. 
In particular, it may lie in a "spin glass" universality class. Understanding the 
nature of color screening and "long range order" in this universal dynamics of- 
fers possibilities for progress in resolving fundamental QCD questions regarding 
properties such as confinement and chiral symmetry. A specific area of progress 
is in the mapping of the CGC degrees of freedom to the traditional language of 
Pomerons and the consequent prospects for understanding soft QCD dynamics. 
In nuclear collisions, CGC dynamics produces "Glasma" field configurations at 
early times: strong longitudinal chromo-electric and chromo-magnetic fields color 
screened on transverse distance scales 1/Q S - These generate long range rapid- 
ity correlations, "sphaleron-like" topological fluctuations characterized by large 
Chern-Simons charge, and instabilities analogous to those seen in QED plasmas. 

The CGC framework is applicable to a wide variety of processes in e+p/A, p+A 
and A+A collisions. It provides an ab initio approach to study thermalization in 
heavy ion collisions and the initial conditions for the evolution of a thermalized 
Quark Gluon Plasma (QGP). The interaction of hard probes with the Glasma is 
little understood and is important for quantifying the transport properties of the 
QGP precursor. A further phenomenological application of the CGC is to the 
physics of ultra high energy cosmic rays (7). 

There are several reviews that discuss various aspects of the CGC and its 
applications in depth (8-10). Our aim here is to provide a broad pedagogical 
overview of the current status of theory and phenomenology for non-experts. 
The next section will focus on the theoretical status, while section 3 discusses the 
CGC in the context of DIS, hadronic and nuclear collisions. Section 4 presents 
the experimental results that probe the dynamics of saturated gluons. We will 
end this review with a brief outlook. 

2 Color Glass Condensate: theoretical status 

We will begin this section with an elementary discussion of QCD bremsstrahlung 
in the Regge-Gribov limit and use this discussion to motivate the BFKL equation, 
how it leads to gluon saturation and describe key features of saturation. The 
theoretical status of the CGC effective theory is then presented. The section 
ends with a discussion of some advanced topics that highlight open issues. 
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2.1 Parton evolution at high energy 



" HMMMM 



Figure 3: Left: elementary bremsstrahlung radiation. Right: high energy scat- 
tering with evolution and recombination. 

The structure of a hadron depends upon the scales resolved by an external 
probe. Quantum fields bound within the hadron can radiate other virtual quanta 
which, by the uncertainty principle, live only for a short period of time. Therefore, 
what looks like an isolated quantum at a given resolution scale reveals a compli- 
cated substructure when probed with finer space-time resolution. The structure 
of the hadron in the IMF is specified by the longitudinal (k + ) and transverse 
(k±) momentum distributions of its quanta and by their correlations-the latter 
become increasingly important with increasing density. 

In perturbative QCD, parton evolution proceeds via bremsstrahlung, which 
favors the emission of soft and collinear gluons. The left part of figure 3 illustrates 
this elementary radiation process. When i < 1, to lowest order in a s , the 
differential probability for this emission is given by 



a s C R d 2 fc ± dx 
d^ rcm - -s- _ - , (1) 



where Cr is the SU(iV c ) Casimir in the color representation of the emitter- N c for 
a gluon and (N^ — 1)/2N C for a quark. This formula demonstrates the collinear 
(k± — > 0) and soft (x — > 0) singularities mentioned above, which produce a 
logarithmic enhancement of gluon emission at small k± and/or x. If the emitter 
at small x were a quark instead a gluon, there would be no small x enhancement; 
while the collinear enhancement is present for either emitter. This asymmetry is 
due to the spin 1 nature of the gluon. 

A high energy scattering such as the one illustrated in the right part of figure 3 
probes the number of gluons with a given x and transverse momenta k± < Q. 
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From eq. (1), the gluon number is 

where the cutoff A QCD has been introduced as a crude way to account for con- 
finement: when confined inside a hadron, a parton has a minimum virtuality of 
0(A 2 ). 

V QCD'' 

In QCD, gluons can radiate softer gluons and thus rapidly multiply as illus- 
trated in figure 3. Each subsequent emission is ot s suppressed however when the 
final value of x is small, these corrections become large. For a cascade with n 
intermediate gluons strongly ordered in x, one obtains 

q „ ^ /'dx^... /"dx: = i / ln iy (3) 

Jx x n J x n x n—l Jx2 X ^ n ' V X J 

When a s ln(l/x) > 1, the correct result for the gluon distribution is obtained by 
summing contributions from all such ladders. The sum exponentiates, modifying 
eq. (2) into 

x dN ° ~ gfgg J_ e "«.y , Y = ln-, (4) 
dxdk 2 ^ tt k\ x 

where w is a number of order unity which is not fixed by this rough estimate. 
The variable Y is known as the rapidity. 

To go beyond this simple power counting argument, one must treat more ac- 
curately the kinematics of the ladder diagrams and include the associated virtual 
corrections. The result is the previously mentioned BFKL equation for the Y- 
evolution of the unintegrated gluon distribution. The solution of this equation, 
which resums perturbative corrections (a s Y) n to all orders, confirms the expo- 
nential increase in eq. (4) , albeit with a /cj_-dependent exponent and modifications 
to the /cj 2 -spectrum of the emitted gluons. 

An important property of the BFKL ladder is coherence in time. Because the 
lifetime of a parton in the IMF, Ax + ~ k + jk\ oc x, the "slow" gluons at the 
lower end of the cascade have a much shorter lifetime than the preceding "fast" 
gluons. Therefore, for the purposes of small x dynamics, fast gluons with x' 3> x 
act as frozen color sources emitting gluons at the scale x. Because these sources 
may overlap in the transverse plane, their color charges add coherently, giving 
rise to a large color charge density. The average color charge density is zero by 
gauge symmetry but fluctuations in the color charge density are nonzero and 
increase rapidly with 1/x. 
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These considerations are at the heart of the CGC reformulation of BFKL evolu- 
tion. However, in contrast to the original BFKL formulation, the CGC formalism 
(11-15) includes non-linear effects which appear when the gluon density becomes 
large. 

The quantity which controls gluon interactions in the IMF is their occupation 
number-the number of gluons of a given color per unit transverse phase-space 
and per unit rapidity, 

n(Ykl bi) - ^ (5) 

where the impact parameter b± is the gluon position in the transverse plane. If 
n<C 1, the system is dilute and gluon interactions are negligible. When n ~ 0(1) 
gluons start overlapping, but their interactions are suppressed by a s <C 1. The 
interaction strength becomes of order one when n ~ 0(l/a s ). It is then that 
non-linear effects become important leading to gluon saturation. 

Gluon occupancy is further amplified if instead of a proton we consider a large 
nucleus with atomic number A > 1. The transverse area scales like A 2 / 3 , and 
the gluon occupation number scales as A 1 / 3 . Thus, for a large nucleus, saturation 
effects become important at larger values of x than for a proton. 

To be more specific, let us discuss a simplified version of the non-linear evo- 
lution equation for the occupation number. Consider an elementary increment 
in rapidity: Y — > Y + dY. Each preexisting gluon in the hadron has a probabil- 
ity a s dY to emit an additional soft gluon - the average increase is dn ~ a s ndY. 
Further, the emission vertex is non-local in k± because the transverse momentum 
of the parent gluon is shared among the two daughter gluons. At high-energies, 
this non-locality is well approximated as a diffusion in the logarithmic variable 
t = ln(k\). Finally, two preexisting gluons can merge and produce a single final 
gluon with rapidity Y + dY. This process is quadratic in a s n and leads to a neg- 
ative term in the evolution equation. Adding up these three effects, one obtains 
the evolution equation 

— ~ ua s n + xa s dt n ~ /3a 2 s n 2 , (6) 

where oj, %■> an d /3 are numbers of order unity. This equation mimics the BFKL 
equation (5) if one drops the term quadratic in a s n in the r.h.s., and mimics its 
non-linear extensions that include saturation if one keeps the quadratic term. 
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2.2 Generic features of gluon saturation 

Although a toy model, eq. (6) captures essential features of saturation 4 . If a s n <C 
1, one can neglect the quadratic term, and eq. (6) predicts an exponential growth 
in Y of the gluon occupation number. But when a s n ~ 1, the negative non- 
linear term turns on and tames the growth. In fact, eq. (6) has a fixed point at 
n = oj/((3a s ) where its r.h.s. vanishes-the evolution stops when this value of n is 
reached, resulting in gluon saturation in the spirit of early works (6). 

Eq. (6) also reveals the emergence of a transverse momentum scale Q S {Y) that 
characterizes saturation. This scale is the k± where gluon occupancy becomes 
of 0(1/ a s ). As a function of transverse momentum, the occupation number 
n(y, is (D(l/a s ) if k± < Q S (Y), and decreases rapidly above Q S (Y) (n oc \jk\ 
for k± » Q S (Y)). The shape of n(Y, k±) as a function of k± is known as the 
saturation front. Note that Q s , the typical gluon transverse momentum at the 
rapidity Y, increases with energy and becomes a semi-hard scale (Q S (Y) » A qCD ) 
at sufficiently high energy. 

Eq. (6) belongs to the generic class of reaction- diffusion processes (19). These 
are processes where an entity can hop to neighboring locations (diffusion term 
Xa s d^n), can split into two identical entities (the term wa s ii), and where two 
entities can merge into a single one (the term f3a^n 2 ). In the limit of large 
occupation numbers, these processes admit the mean field description of eq. (6), 
which in the context of statistical physics is known as the Fisher-Kolmogorov- 
Petrovsky-Piscounov (FKPP) equation (19). In QCD, the closest equation of this 
type is the Balitsky-Kovchegov (BK) equation (13). The correspondence between 
the BK and FKPP equations, originally noticed in (20), clarifies the properties of 
saturation fronts in QCD in analogy with known properties of reaction-diffusion 
processes. 

A crucial property is the emergence of traveling waves. The saturation front 
generated by this equation propagates without distortion at constant speed; one 
has n(Y, t) = n(t — \ S Y) with X s a constant. This property has been verified 

4 Our discussion oversimplifies the mechanism for gluon saturation. In the saturation regime, 
gluons form configurations that screen their color charge over transverse scales ~ 1/Q S (16,17). 
Thus, gluons with momentum k± Q s are emitted from a quasi-neutral patch of color sources, 
and their occupation number grows only linearly in Y (16,18) - much slower than the exponential 
growth in the region k± 3> Q s . 
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in numerical studies of the BK equation (21, 22) and in analytic studies of the 
BFKL equation in the presence of a saturation boundary (23,24). It provides a 
natural explanation of the geometric scaling phenomenon observed in the HERA 
data (25,26) (see section 3.1). In QCD, a front moving with constant speed A s 
is equivalent to the saturation momentum increasing exponentially with Y, 

Q 2 s (Y)~Q 2 e x ° Y with A s «4.9a s , (7) 

where Qq is some non-perturbative initial scale. For a large nucleus, Qq scales 
like A 1 / 3 as does Q 2 (Y) for any Y. This form of the saturation momentum is 
modified to Q 2 = Qq e^ x ^ Y+Y °^ when the running of the strong coupling is taken 
into account; see (24) for a detailed study of higher order effects on the energy 
dependence of Q s . 



2.3 The Color Glass Condensate 

The CGC is an effective field theory (EFT) based on the separation of the degrees 
of freedom into fast frozen color sources and slow dynamical color fields (11). A 
renormalization group equation -the JIMWLK equation (14, 15)- ensures the 
independence of physical quantities with respect to the cutoff that separates the 
two kinds of degrees of freedom. 

The fast gluons with longitudinal momentum k + > A + are frozen by Lorentz 
time dilation in configurations specified by a color current Jtt = p a , where 
p a (x~,x±) is the corresponding color charge density. On the other hand, slow 
gluons with k + < A + are described by the usual gauge fields of QCD. Because 
of the hierarchy in k + between these two types of degrees of freedom, they are 
coupled eikonaly by a term J^A^. The fast gluons thus act as sources for the fields 
that represent the slow gluons. Although it is frozen for the duration of a given 
collision, the color source density p a varies randomly event by event. The CGC 
provides a gauge invariant distribution W A + [p] , which gives the probability of a 
configuration p. This functional encodes all the correlations of the color charge 
density at the cutoff scale A + separating the fast and slow degrees of freedom. 
Given this statistical distribution, the expectation value of an operator at the 
scale A + is given by 

(0) A+ ^ J [Dp] W A+ [p] O [p] , (8) 

where 0[p\ is the expectation value of the operator for a particular configuration 
p of the color sources. 
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The power counting of the CGC EFT is such that in the saturated regime the 
sources p are of order g~ l . Attaching an additional source to a given Feynman 
graph does not alter its order in g; the vertex where this new source attaches to the 
graph is compensated by the g~ l of the source. Thus, computing an observable 
at a certain order in g 2 requires the resummation of all the contributions obtained 
by adding extra sources to the relevant graphs. The leading order in g 2 is given by 
a sum of tree diagrams, which can be expressed in terms of classical solutions of 
the Yang-Mills equations. Moreover, for inclusive observables (27), these classical 
fields obey a simple boundary condition: they vanish when t — > — oo. 

Next-to-leading order (NLO) computations in the CGC EFT involve a sum 
of one-loop diagrams embedded in the above classical field. To prevent double 
counting, momenta in loops are required to be below the cutoff A + . This leads 
to a logarithmic dependence in A + of these loop corrections. These logarithms 
are large if A + is well above the typical longitudinal momentum scale of the 
observable considered, and must be resummed. 

For gluon correlations inside the the hadron wavefunction and also for suffi- 
ciently inclusive observables in a collision, the leading logarithms are universal 
and can be absorbed into a redefinition of the distribution [p] of the hard 
sources. The evolution of [p] with A + is governed by the functional JIMWLK 
equation 

where % is known as the JIMWLK Hamiltonian. This operator contains up to 
two derivatives d/dp, and arbitrary powers in p. Its explicit expression can be 
found in refs. (8,14,15). The derivation of the JIMWLK equation will be sketched 
in the section 3.1. 

Numerical studies of JIMWLK evolution were performed in (22,28). An an- 
alytic, albeit formal, solution to the JIMWLK equation was constructed in (29) 
in the form of a path integral. Alternatively, the evolution can can be expressed 
as an infinite hierarchy of coupled non-linear equations for n-point Wilson line 
correlators-often called the Balitsky hierarchy (12). In this framework, the BK 
equation is a mean field approximation of the JIMWLK evolution, valid in the 
limit of a large number of colors N c — > oo. Numerical studies of the JIMWLK 
equation (22, 28) have found only small differences with the BK equation. 

Let us finally comment on the initial condition for the JIMWLK equation 
which is also important in understanding its derivation. The evolution should 
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start at some cutoff value in the longitudinal momentum scale Aq" at which the 
saturation scale is already a (semi)hard scale, say Q s o > 1 GeV, for perturbation 
theory to be applicable. The gluon distribution at the starting scale is in general 
non-perturbative and requires a model. A physically motivated model for the 
gluon distribution in a large nucleus is the McLerran-Venugopalan model (11). 
In a large nucleus, there is a window in rapidity where evolution effects are not 
large but x is still sufficiently small for a probe not to resolve the longitudinal 
extent of the nucleus. In this case, the probe "sees" a large number of color 
charges, proportional to A 1 / 3 . These charges add up to form a higher dimensional 
representation of the gauge group, and can therefore be treated as classical color 
distributions (11,30). Further, the color charge distribution W A + [p] is a Gaussian 
distribution 5 in p. The variance of this distribution -the color charge squared 
per unit area- is proportional to A 1 / 3 and provides a semi-hard scale that makes 
weak coupling computations feasible. In addition to its role in motivating the 
EFT and serving as the initial condition in JIMWLK evolution, the MV model 
allows for direct phenomenological studies in p+A and A+A collisions in regimes 
where the values of x are not so small as to require evolution. 

The dynamics of small x gluons in QCD may be universal in more than one 
sense (33). A weak form of this universality is that their dynamics in both hadrons 
and nuclei is controlled only by the saturation scale with its particular dependence 
on energy and nuclear size. A stronger form of the universality is noticed in 
particular for the solution of the BK equation with running coupling effects; the 
saturation scale, for both hadrons and nuclei, at fixed impact parameter, becomes 
the same asymptotically with increasing energy (34). The strongest form of the 
universality is that the RG flows in the saturation regime have a fixed point 
corresponding to universal "critical" exponents describing the behavior of multi- 
parton correlation functions. As discussed further below, the RG equations for 
high energy QCD lie in a wide class of reaction-diffusion processes which have 
universal properties remarkably close to those of spin glasses (35). 

5 There is a additional term, corresponding to the cubic Casimir; which is parametrically 

suppressed for large nuclei (31). This term generates Odderon excitations in the JIMWLK/BK 
evolution (32). 
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2.4 Advanced Theory topics 

Reaction-diffusion processes exhibit an extreme sensitivity to particle number 
fluctuations (36-38), generated by gluon splittings, which produce correlations 
among pairs of gluons (39) . This effect is of higher order in a s , and is linear in 
n since it results from the splitting of a single gluon. Conversely, producing two 
gluons without this splitting leads to a term that has one less power of a s , but 
of order n 2 . Thus, the splitting contribution is important in the dilute regime 
where n < a s . Since, as mentioned previously, the dynamics of the saturation 
front is driven by the BFKL growth of its dilute tail, these fluctuations play an 
important role in the 2-gluon density (nn) at high energy. 

As manifest in eq. (6), this 2-gluon density enters in the non-linear term lead- 
ing to saturation of the single gluon density. Thus, gluon number fluctuations 
in the dilute regime can strongly influence the approach towards saturation. For 
instance, as argued in (36) for generic reaction-diffusion processes, and indepen- 
dently in the QCD context(37), these fluctuations reduce the (average) speed of 
the saturation front. Besides making the value of Q s a fluctuating quantity, they 
tend to wash out the geometric scaling property of the individual fronts (38,39). 
Both effects are quantitatively important, as shown by explicit numerical simu- 
lations within various reaction-diffusion models (including those inspired by the 
QCD dynamics at high energy (40)). 

There is presently no general theory that includes BFKL ladders, saturation 
and fluctuations. (In terms of Feynman graphs, this corresponds to resumming 
"Pomeron loops" diagrams to all orders (39).) Attempts to construct such a 
theory (39,41,42) have led to incomplete formalisms that are difficult to exploit 
in phenomenological applications. Fortunately, the effect of these fluctuations is 
considerably reduced by the running of the strong coupling (43), which tends to 
postpone their importance to unrealistically large rapidities. Thus, for practical 
applications at least up to LHC energies, a sufficient theory for the approach 
towards saturation is the leading-order mean-field evolution extended with run- 
ning coupling corrections. This theory has developed significantly as the running- 
coupling version of the BK equation has been constructed in (44) and successfully 
applied to studies of the phenomenology at HERA (45). 
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3 Collisions in the CGC framework 

The CGC is an effective theory for the wavefunction of a high-energy hadron or 
nucleus. In this section, we apply it, with particular emphasis on factorization, 
to deeply inelastic scattering and hadronic collisions. In A+A collisions, the 
formation of the Glasma and its key features are emphasized. 

3.1 The CGC and DIS at small x 

At small x, corresponding to large Ioffe times (46), DIS is characterized by the 
fluctuation of the virtual photon into a quark-antiquark pair which then scat- 
ters off the hadronic or nuclear target. The inclusive DIS cross-section can be 
expressed as (47) 



where ip(z,r±) is the qq component of the wave-function of the virtual photon 
(known from QED) and CdipoieC^ is the QCD "dipole" cross-section for the 
quark-antiquark pair to scatter off the target. This process is shown in figure 4, 
where we have assumed that the target moves in the — z direction. In the leading 
order (LO) CGC description of DIS, the target is described, as illustrated in 
figure 5, as static sources with k~ > Aq . The field modes do not contribute at 
this order. 



Figure 4: Left: leading Order (LO) contribution to DIS off the CGC. Right: NLO 
contribution. 

Employing the optical theorem, a,n P o\e(x, r±) can be expressed in terms of the 
forward scattering amplitude T(x±, y±) of the qq pair at LO as 




(10) 





(11) 




(12) 
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with U(x±) a Wilson line representing the interaction between a quark and the 
color fields of the target, defined to be 

rl/xP- 

U(x±) = T expig / dz + A~(z + , x±) . (13) 



In this formula, A" is the minus component of the gauge field generated (in 
Lorenz gauge) by the sources of the target; it is obtained by solving classical 
Yang-Mills equations with these sources. The upper bound xP~ (where P~ 
is the target longitudinal momentum and x the kinematic variable defined in 
figure 1) indicates that source modes with k~ < xP~ do not contribute to this 
scattering amplitude. Thus if the cutoff Aq of the CGC EFT is lower than xP~ , 
T LO is independent of Aq~. 

However, when Aq is larger than xP~ , the dipole cross-section is in fact inde- 
pendent of x (since the CGC EFT does not have source modes near the upper 
bound xP~) and depends on the unphysical parameter An. As we shall see now, 
this is related to the fact that eq. (11) is incomplete and receives large correc- 
tions from higher order diagrams. Consider now the NLO contributions (one of 
them is shown in the right panel in figure 4) with gauge field modes in the slice 
Aj~ < k~ < Aq (see figure 5). An explicit computation of the contribution of 

fields »< — sources — » 



AS 



+ 



fields 



H — I— 

A] Ag 



5T T 

NLO U 



Figure 5: Top: sources and fields in the CGC effective theory. Bottom: NLO 
correction from a layer of field modes just below the cutoff. 



field modes in this slice gives 

*T NLO ( aJ± ,yJ=ln^) UT LO (x ± ,y ± ), (14) 

where % is the JIMWLK Hamiltonian. All dependence on the cutoff scales is 
in the logarithmic prefactor alone. This Hamiltonian has two derivatives with 
respect to the classical field A ~ 0(1/ g); HT LO is of order a s T LO and there- 
fore clearly an NLO contribution. However, if the new scale A^~ is such that 
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ct s ln(A^"/A^) ~ 1, this NLO term becomes comparable in magnitude to the LO 
contribution. Averaging the sum of the LO and NLO contributions over the 
distribution of sources at the scale Aq , one obtains 



J [Dp] W A -[p] (T LO +5T NLO ) = J [Dp] W A -[p] T L 



(15) 



where W A - = (1 + ln(A /A l )H)W A -. We have shown here that the NLO 
correction from quantum modes in the slice A^ < k~ < Aq can be absorbed in 
the LO term, provided we now use a CGC effective theory at A^~ with the modified 
distribution of sources shown in eq. (15). In differential form, the evolution 
equation of the source distribution, 

jj^zyWA- ~«W A -, (16) 

is the JIMWLK equation. 

Repeating this elementary step, one progressively resums quantum fluctuations 
down to the scale k~ ~ xP~ . Thanks to eq. (15), the result of this resummation 
for the dipole cross-section is formally identical to eq. (11), except that the source 
distribution is W x p- instead of W A - . Note that if one further lowers the cutoff 
below xP~, the dipole cross-section remains unchanged. 



3.2 The CGC in p+A collisions 





Figure 6: Left: sketch of a proton-nucleus collision. Right: example of leading 
order contribution in a nucleus-nucleus collision. 



Collisions between a dilute hadron projectile and a dense hadron target can be 
studied semi-analytically in the CGC framework. The archetype of such collisions 
is a proton-nucleus collision. However, the dilute-dense treatment also applies to 
proton-proton collisions for measurements at forward rapidities where the wave- 
function of one of the projectiles is probed at large x and that of the other at 
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small x. These asymmetrical collisions can be treated using conventional parton 
distributions for the proton and the CGC for the nucleus. 

The simplest quantity to compute in this context is the single particle inclusive 
spectrum. In the CGC framework, this process is due to the scattering of a 
parton from the proton off the color field of the target nucleus, as illustrated in 
figure 6 (left). When the parton is a quark, the amplitude for the scattering is 
proportional to the Fourier transform of a Wilson line, 

M LO oc J d 2 x ± e ip ^ x± U(x±) . (17) 

Squaring this amplitude, summing over the color of the quark in the final state 
and averaging over the color of the incoming quark, we get 

\M\l Q oc J d*b d"r ± e*^ T LO (b + ^, b - ^) , (18) 

where T LO is the dipole scattering amplitude (see eq. (12)) that already appeared 
in our discussion of DIS. Because the same quantity appears here, the treatment 
of the NLO corrections we described in the DIS case is similarly applicable; one 
integrates out softer modes by lowering the cutoff of the CGC EFT letting W\- [p] 
evolve according to the JIMWLK equation. The scale to which one evolves 
the cutoff is A~ = xP" with x = (p±/y/s) exp(— y) where p± is the transverse 
momentum of the scattered parton and y its rapidity. Therefore, in the CGC 
framework, the cross-section for this process is simply the Fourier transform of 
the dipole cross-section (T^i po \ e (x, r±). In p+A collisions, final states containing 
a photon or a lepton pair can be similarly expressed in terms of the same Fourier 
transform (50). 

An additional contribution to the single inclusive particle spectrum in a p+A 
collision is due to an incoming gluon instead of a quark. The treatment is nearly 
identical to the incoming quark case, except that in eq. (17) one must replace 
the Wilson line in the fundamental representation by a Wilson line in the adjoint 
representation. 

Similar calculations can be performed for processes with more complicated final 
states, such as the production of a quark-antiquark pair. Although this observable 
has a more complicated expression (containing terms that are the product of four 
Wilson lines) its NLO corrections still comply with eq. (14), which ensures their 
factorization into the distribution of sources WX-[p]. The crucial ingredient for 
factorization to work is to consider an observable that is sufficiently inclusive to 
allow the corresponding final state to be accompanied by an arbitrary number 
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of gluons. Any restriction on associated gluon radiation will not permit NLO 
corrections to factorize simply in the distribution of sources. 

A much more limited but widely used form of factorization is k± factoriza- 
tion (51) in terms of k± dependent unintegrated quark and gluon distributions 
of the projectile and target. Within the CGC framework, these results can be re- 
produced for gluon (52,53) and heavy quark (54) distributions at large transverse 
momenta kj_ > Q s ; however, at smaller transverse momenta, k± factorization is 
broken even at leading order (55-57). 



3.3 Shattering CGCs in A+A collisions 

Collisions between two nuclei ("dense-dense" scattering) are complicated to han- 
dle on the surface. However, in the CGC framework, because the wave functions 
of the two nuclei are saturated, the collision can be treated as the collision of 
classical fields. This insight significantly simplifies the treatment of A+A scat- 
tering. The classical fields are coupled to fast partons of each nucleus respectively 
described by the external current = p\ + pi. The source densities of 
fast partons p± t 2 are both parametrically of order 1/g, which implies that graphs 
involving multiple sources from both projectiles must be resummed. (See the 
right panel of figure 6 for an illustration.) 

At leading order, inclusive observables 6 depends on the retarded classical color 
field A^, which solves the Yang-Mills equations [P„, J 7 ^] = J u with the boundary 
condition linL E o^_ 00 A^ 1 = 0. Among the observables to which this result applies 
is the expectation value of the energy-momentum tensor at early times after the 
collision. At leading order, 

T Z = -^f ~ T^F"* , (19) 

where F^ u is the field strength of the classical field A^. 

Although A+A collisions are more complicated than e+A or p+A collisions, 
one can still factorize the leading higher order corrections into the evolved dis- 
tributions Wa-[pi] and Wa+[/»2]- At the heart of this factorization is a gener- 
alization of eq. (14) to the case where the two projectiles are described in the 
CGC framework (58). When one integrates out the field modes in the slices 

6 Exclusive observables may also be expressed in terms of solutions of the same Yang-Mills 
equations, but with more complicated boundary conditions than for inclusive observables. 
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Af < < A^, the leading correction to the energy momentum tensor is 



ST**" 



NLO 



h '£) W ' +ta ($'* 



T% , (20) 



where Hi,2 are the JIMWLK Hamiltonians of the two nuclei respectively. What 
is crucial here is the absence of mixing between the coefficients of the loga- 
rithms of the two projectiles; they depend only on p\^ respectively and not on the 
sources of the other projectile. Although the proof of this expression is somewhat 
involved, the absence of mixing is deeply rooted in causality. The central point 
is that because the duration of the collision (which scales as the inverse of the 
energy) is so brief, soft radiation must occur before the two nuclei are in causal 
contact. Thus logarithms associated with this radiation must have coefficients 
that do not mix the sources of the two projectiles. 

Following the same procedure for eq. (20), as for the e+A and p+A cases, one 
obtains for the energy-momentum tensor in an A+A collision the expression 

<^> LLog = J [D Pl D P2 ] Wl [ Pl ] W 2 [p 2 ] T% . (21) 

This result can be generalized to multi-point correlations of the energy-momen- 
tum tensor, 

(T^( Xl )...T^(x n )) LLog = j[Dp l Dp 2 ]W 1 [ Pl ]W 2 [p 2 ] 

xT^( Xl )---T^(x n ). (22) 

In this expression, all the correlations between the energy-momentum tensor at 
different points are from the distributions ^1,2^1,2]- Thus, the leading corre- 
lations are already built into the wavefunctions of the projectiles prior to the 
collision. 

The expressions in eqs. (21) and (22) are valid for proper times r ~ 1/Q S after 
the heavy ion collision. The energy-momentum tensor, for each configuration of 
sources pip is determined by solving classical Yang-Mills equations to compute 
the gauge fields A c ^[p\, p 2 ] in the forward light cone with initial conditions deter- 
mined by the classical CGC fields of each of the nuclei at r = (55,59-62). The 
corresponding non-equilibrium matter, with high occupation numbers ~ l/ct s 
is called the Glasma (63). The Glasma fields at early times are longitudinal 
chromo-electric and chromo- magnetic fields that are screened at distances 1/Q S 
in the transverse plane of the collision. As a consequence, the matter produced 
can be visualized (see figure 7) as comprising R\Q 2 S color flux tubes of size l/Q s , 
each producing l/a s particles per unit rapidity. The flux tube picture is sup- 
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Figure 7: Left: Gauge field configurations in the form of "flux tubes" of longi- 
tudinal chromo-electric and chromo-magnetic fields screened on transverse scales 
1 IQ S . Right: Model comparison (64 ) to long range rapidity correlations measured 
by the PHOBOS collaboration (65). 

ported by non-perturbative numerical solutions of the classical Yang- Mills equa- 
tions (66). The "Glasma flux tubes" generate n-particle long range rapidity 
correlations (58,64,67,68). These distributions are negative binomial distribu- 
tions (69). They also carry topological charge (70); the resulting dynamical topo- 
logical "sphaleron" transitions may result in observable metastable CP-violating 
domains (71). 

The evolution of the Glasma into a thermalized Quark Gluon Plasma (QGP) 
is not understood. An important ingredient is the role of instabilities (72). At 
early times, these arise at NLO from terms that break the boost invariance of the 
LO term (10,58). The modification to the evolution of the Glasma is obtained by 
solving 3+1-D Yang-Mills equations (73) for the (now) rapidity dependent gauge 
fields convolved with a distribution giving the spectrum of fluctuations (74). 
While these effects may isotropize the system, early thermalization may also 
require collisions whose role still needs to be clarified (75). 

4 Phenomenological applications of the CGC 

In this section, we will discuss the applications of the theoretical formalism out- 
lined in the previous section to analyze and predict a wide range of phenomena 
ranging from DIS in e+p and e+A collisions to the scattering of hadronic pro- 
jectiles ranging from p+p to p+A to A+A collisions. A unifying ingredient in 
many of the applications is the dipole cross-section defined in eq. (11), albeit, 
as apparent in the treatment of A+A collisions, the fundamental ingredient is 
really the density matrix Wy[p]. Because the JIMWLK equation (eq. (16)) for 
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this quantity is time consuming to solve 7 , many of the applications are in the 
context of models of the dipole cross-section which incorporate key features of 
saturation. These models provide an economical description of a wide range of 
data with only a few parameters. A good compromise between the full JIMWLK 
dynamics and models of the dipole cross-section is the BK equation, which is a 
large N c realization of JIMWLK dynamics. With the recent availability of the 
NLO BK equation, global analyzes of data are in order. Much of our discussion 
below is in the context of dipole models; improvements a la BK are highlighted 
wherever available. 



4.1 DIS in e+p and e+A collisions 

A remarkable observation (25) is that HERA data (25,26) on the inclusive virtual 
photon-proton cross section for x < 0.01 scale as a function of the ratio Q 2 /Q 2 (x); 
see the left part of figure 8. This scaling is violated for larger values of x. Here 
Q 2 S = Qoi x o/%) X is the saturation scale with Qq = 1 GeV 2 , xq = 3 • and 
A ~ 0.3. This scaling is referred to as "geometrical scaling", because the survival 
probability of the color dipole that the virtual photon fluctuates into is close to 
unity or zero respectively depending on whether the ratio of the saturation radius 
(~ 1/Qs) to the size of the dipole (of size ~ 1/Q) is large or small. Recall that 
the saturation radius denotes the typical size of regions with strong color fields. 
Geometric scaling has also been observed in inclusive diffraction, exclusive vector 
meson production and deeply virtual Compton scattering data at HERA (78). In 
detail, the data also show violations of geometric scaling, which can be interpreted 
as consequences of BFKL diffusion (79), non-zero quark masses (80) and possibly 
DGLAP evolution as well (81). Note that the best scaling is obtained with a 
saturation scale that behaves like Q 2 s {x) oc x -0,3 , a slower x-dependence than 
predicted by the LO BK equation. This discrepancy is resolved by a resummed 
NLO computation of the saturation exponent (24) which indeed gives 0.3. 

While geometrical scaling is very suggestive of the presence of semi-hard dy- 
namical scales in the proton, it is not conclusive in and of itself (82); more detailed 
comparisons to the data are essential. Despite their simplicity, saturation mod- 
els (45,79-81,83-86) provide remarkably good descriptions of HERA data at small 
x < 0.01. The free parameters are fixed from fits to the total cross-section data 

7 Other unknowns include higher order corrections, initial conditions at low energy and 

impact-parameter dependence of distributions. 
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Figure 8: Geometrical scaling and limiting fragmentation. Left: <r 7 * p data at 
HERA for x < 0.01 and all Q 2 up to 450 GeV 2 ; r is the scaling variable, r = 
Q 2 /Q 2 (x) (25,26). Right: particle multiplicities for several collision energies at 
RHIC (76), compared to the computation of (77). 

alone; once these are fixed, the models predict a large variety of results, includ- 
ing the longitudinal (F L ), diffractive (F^), an< ^ charm (Ffi) structure functions, 
the virtual photon production of vector mesons (p, J/ip), and the deeply virtual 
Compton scattering (DVCS). The most recent analysis of inclusive data (45) of 
DIS in e+p collisions is quite sophisticated; the energy dependence is given by 
the running-coupling BK equation, and the free parameters refer solely to the 
initial conditions and to the proton transverse area. 

The phenomenon of hard diffraction in DIS is particularly sensitive to satura- 
tion. The simplest diffractive processes are events in which the proton remains 
intact and a large gap in rapidity with no particles extends between the rapidity 
of the proton and that of the fragmentation products of the virtual photon. For 
small invariant masses, this process corresponds to elastic scattering of the qq 
dipole off the target. Its cross-section is evaluated as 8 

l 

^lf| _ = lo^/^/ d2 rJVK^±)| 2 ^dipole(^±) • (23) 

The dipole cross-section in this expression, for small dipoles r 2 <C 1/Q 2 , is a 
color singlet combination of two gluons that can be interpreted as Pomeron ex- 
change (87); for larger dipoles r 2 > 1/Q 2 , the color singlet exchange does not 
have this simple interpretation. The t distribution has the form dadis/dt = 
8 Forward diffraction (corresponding to t = 0, where t — (P — P') 2 is the momentum transfer 

squared between the incoming and outgoing proton) can be compared directly to inclusive DIS 

within a dipole model since it depends only on the dipole cross-section. 
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Figure 9: Left: ZEUS data for the ratio a dm /ctot together with the corresponding 
prediction of the saturation model in (84). Right: DVCS measurement at HERA, 
and comparison with (81). 

exp(—B D \t\)dadis/dt\t=o, where B D is the transverse area of the interaction re- 
gion in the proton and is closely related to the transverse gluon radius in the 
proton estimated to be 0.61 ± 0.04 fm (88). From this form of the diffractive 
cross-section and eq. (23), the total diffractive cross-section is 




We used here the color transparency approximation <7di po ie oc r 2 Q 2 s {x) for the 
dipole cross-section. Unlike the inclusive cross-section in eq. (10) which is domi- 
nated by small dipole sizes ~ l/Q, the integrand of the diffractive cross-section 
is dominated by larger size dipoles of size ~ 1/Q S - Comparing eq. (24) with 
eq. (10), we deduce that in the saturation framework the ratio Odift/otot * s a P~ 
proximately constant as a function of energy. As shown in figure 9 (left), the 
HERA data support this qualitative observation and are in quantitative agree- 
ment with the detailed saturation model of (84). We note further that excellent 
fits with x 2 ~ 1 are obtained in this saturation framework for exclusive vector 
meson production and deeply virtual Compton scattering (81) in addition to in- 
clusive diffraction (89) (see figure 9, right). These exclusive processes provide 
detailed information about the impact parameter dependence of the dipole cross- 
section (90) and may even provide unique information about the partonic nature 
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of short range nuclear forces (88). 

In contrast, if unitarization were due to soft physics (the prevailing viewpoint 
before the advent of saturation (91)), the natural cutoff for the integral would be 
1 / A QCD . Diffraction would be non-perturbative even for hard Q and the energy 
dependence of the diffractive cross-section would be the square of the inclusive 
cross-section, in disagreement with data. 

Fits based on geometrical scaling of the e+A fixed target data with the func- 
tional form Ql oc A 5 (with 5 naively 1/3) give values of 5 with range from 1/4 
to 4/9 (92). However, the dipole formalism which describes e+p data so success- 
fully can be straightforwardly generalized to nuclei to construct the corresponding 
nuclear dipole cross-section (85). When compared to results for the x,Q 2 depen- 
dence of inclusive e+A data, one finds effectively that Q* >A (x) = Q^(x)^ 1 / 3 (93). 
Nuclear diffractive distributions can be computed in saturation models (89,94); 
predictions have been made as well for semi- inclusive hadron production (95), ex- 
clusive vector meson production (96) and nuclear DVCS (97) in this framework. 

In saturation models of nuclei, the small x distributions in a nucleon are con- 
volved with nuclear geometry to give the nuclear distributions. This process 
however does not commute with the RG evolution in x of nuclear distributions 
determined at some initial scale. First computations of e+A inclusive distri- 
butions in the NLO BK framework have been performed and good agreement 
obtained for existing fixed target data (64). 

4.2 Particle multiplicities in d+A and A+A collisions 

The CGC EFT is most reliable when at least one of the projectiles is dense in the 
sense discussed in the previous section. At RHIC, the world's first deuteron+hea- 
vy nucleus (d+A) and A+A collider, many features of the CGC are being tested. 
These include bulk features such as the rapidity and centrality dependence of 
particle multiplicities in d+A and A+A collisions, limiting fragmentation, particle 
spectra and correlations, and even possibly more exotic features such as long 
range rapidity correlations ( "the ridge" ) in A+A collisions and local CP violation 
arising from topological fluctuations in A+A collisions. In this sub-section, we 
will focus on particle multiplicities and discuss the other features in subsequent 
sub-sections. 

"Limiting fragmentation" is the well known property of the strong interactions 
that the rapidity distribution in the fragmentation region becomes independent 
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of the collision energy. When one increases the beam energy (see figure 8, right), 
dN/drj' is the same at large rf for all energies, where rf = r] — r/bcam- I n the 
fragmentation region, large x\ modes are probed in one hadron or nucleus and 
small X2 modes in the other. At small X2, if gluonic matter is saturated, parton 
distributions have a very weak dependence on or equivalently on t\ + 7/beam! 
cross-sections will only depend on x\ or rj — r/bcam- Deviations from limiting 
fragmentation at high energies are especially interesting because of a significant 
window in phase space where the RG evolution in X2 to (or away) from the uni- 
versal "black disk" of saturated gluon matter can be explored (77,98). These 
constraints lead to specific predictions for the rapidity dependence of the multi- 
plicity at the LHC in the BK RG framework (77,99). 

The CGC can only predict the distribution of initial gluons, at a proper time 
r ~ Qj 1 - In p/d+A collisions this is not a significant limitation because only 
a few gluons are produced in the final state. The situation is vastly different 
in A+A collisions where the particle multiplicity is significantly larger and the 
system evolves from the non-equilibrium Glasma to a QGP, the latter evolving 
subsequently as a fluid with low viscosity. For p+A collisions, the problem is 
solvable analytically (52,53,100). For A+A collisions, two kinds of calculations 
have been performed: 

(a) Exact numerical solutions of the Yang-Mills equations (55,59,61). Quantum 
evolution effects are not included systematically and the ensemble of color charges 
is assumed to be Gaussian (MV model), which is reasonable at RHIC energies of 
x ~ 10~ 2 at RHIC. The inclusion of quantum effects from the wavefunctions (58, 
69) is under control, but is still an unsolved problem for those present in the final 
state evolution (73,74). 

(b) Approximate analytical calculations of the initial gluon spectrum (101). 
These calculations assume fcj_-factorization although it is violated for momenta 
P± < Qs- In this model, saturation effects are introduced via the unintegrated 
gluon distribution of the nuclei with the rapidity dependence of the gluon spec- 
trum governed by that of the saturation scale Q 2 s {x) ~ x -0 3 . In the CGC frame- 
work, dN/dn = Q 2 s R 2 A /a s (Q s )] a unique feature is that the coupling runs as a 
function of Q s (102), an observation which is in agreement with the centrality 
dependence of RHIC data (103). Similar analyzes of multiplicity distributions 
were extended to the case of d+A collisions (104), with results in fair agreement 
with RHIC data. 
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4.3 Inclusive spectra in d+A and A+A collisions 

4.3.1 dependence One does not expect significant final state inter- 
actions with on-shell partons in p + A collisions. Moreover, the inclusive gluon 
spectrum depends only on the Wilson line correlator (U (0)W(x±)) (53,100,105), 
that can be determined from the dipole cross-section used in DIS (106). The 
hadron spectrum in d+A collisions has been evaluated in this approach in (107), 
with results in good agreement with RHIC data. Very recently, good agreement 
of single inclusive distributions in d+A collisions were obtained in the NLO RG 
approach (108) consistent with its application in e+A collisions (64). 

In A+A collisions, the strong final state interactions likely alter the momentum 
distribution to be more isotropic at low momenta; the hard tail is modified by 
parton energy loss. While these effects may in part be included in the Glasma, 
the additional contributions at later stages are important. 

4.3.2 Nuclear modification ratios To quantify nuclear effects, p± spec- 
tra in A+A and d+A can be compared to those in p+p collisions by the ratios 
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where N co \\ ~ A 4 / 3 (resp. A) is the number of binary collisions in A+A (resp. 
d+A) collisions. At high pj_, one expects the ratios to scale with N co \\. In A+A 
collisions at RHIC, because no suppression is seen in d+A collisions at central 
rapidities, the large observed suppression in Raa is rightly interpreted as a final 
state effect, with a strong candidate being energy loss induced by the dense 
medium (109). 

However, at small x, RG evolution a la BK predicts (110) that this ratio 
should be suppressed by saturation effects in the wavefunction of the nucleus. 
This can be studied in the clean environment of d+A collisions where small x is 
probed measuring the ratio at a positive rapidity in the direction of the deuterium 
nucleus. The ratio R^a was measured by the BRAHMS (111) collaboration up to 
rapidities of rj = 3.2 and by the STAR collaboration (112) at r] = 4. At these large 
rj, the value of RdA is significantly below unity, consistent with predicted trend. In 
addition, as anticipated, the suppression is greater for more central collisions. The 
onset of this forward suppression was also studied semi-analytically in (100,113), 
and more quantitative calculations have been performed in (107,114,115), using 
models of the dipole cross-section that have a realistic x-dependence and most 
recently in the NLO BK framework (108). 
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4.4 Two hadron correlations in d+A and A+A collisions 
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Figure 10: Left: causal relations between two particles separated in rapidity. 
Right: 2-hadron correlation measured by STAR as a function of An and A<j>. 

Two hadron correlations 9 are more sensitive in distinguishing between model 
predictions than the single inclusive results for which alternative explanations of 
the data are feasible. Very recently, striking preliminary results on these cor- 
relations in d+Au collisions have been presented by the STAR and PHENIX 
collaborations. In the STAR measurements, back-to-back correlations of pairs 
have been studied where i) one particle in the pair is forward in rapidity and the 
other at central rapidity and ii) both particles in the pair are at forward rapidi- 
ties. The "forward-central" (FC) pairs probe x ~ 10 -2 in the gold nucleus while 
the "forward- forward" (FF) pairs probe x ~ 10 -3 in the gold nucleus. A clear 
broadening of the backward hadron peak is seen in the transition from FC to FF, 
as well as with increasing centrality in the FF events-in the latter case, the distri- 
bution is so broad that no peak is visible! This particular effect was predicted in 
the CGC (117) and is in quantitative agreement the STAR results (118). These 
results are also consistent with measurements by the PHENIX collaboration on 
pair distributions in d+Au collisions. Related discussions can be found in (119). 

While two particle correlations in A+A collisions are typically much altered 
by final state interactions, there is an important exception for particles widely 
separated in rapidity Causality implies that these correlations are created at 
very early times, as illustrated in the left panel of figure 10. A simple estimate 
for ultra-relativistic particles whose space-time and momentum rapidities are 
strongly correlated gives 



''"max Tfreezeout 6 



(26) 
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Hadron-photon correlations have also been studied in (116). 
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for the latest time at which these particles could have been correlated. For a 
freeze-out time Tf reezC out ~ 10 fm/c, and rapidity separations Arj > 4, one sees 
that these correlations must have been generated before 1 fm/c. Thus long range 
rapidity correlations are a "chronometer" of the evolution of the very strong 
Glasma color fields produced early in the heavy ion collision. 

Striking long range rapidity correlations were seen in A+A collisions in events 
with prominent jet like structures, where the spectrum of associated particles is 
observed to be collimated in the azimuthal separation Aip relative to the jet and 
shows a nearly constant amplitude in the strength of the pseudo-rapidity correla- 
tion At] up to At] ~ 1.5 (120). These events were coined "ridge" events following 
from the visual appearance of these structures as an extended mountain ridge in 
the A-q-Aip plane associated with a narrow jet peak. This collimated correlation 
persists up to Ar] ~ 4 (65,121). An important feature of ridge correlations is that 
the above described structure is seen in two particle correlations without a jet 
trigger (see figure 10) and persists without significant modification for the trig- 
gered events (122,123). These events include all hadrons with momenta p± > 150 
MeV. In such events, a sharp rise in the amplitude of the ridge is seen (122, 124) 
in going from peripheral to central collisions. 

The ridge structures seen in two particle and three particle (125) A+A cor- 
relations can be explained as resulting from the transverse radial flow of the 
Glasma flux tubes we discussed previously in section 3.3. The flux tubes are 
responsible for the long range rapidity correlation in the ridge; the angular col- 
limation occurs because particles produced isotropically in a given flux tube are 
collimated by the radial outward "Hubble" hydrodynamic flow of the flux tubes. 
Ideas on the angular collimation of particle distributions by flow were discussed 
in refs. (126). When combined with the long range rapidity correlations provided 
by the Glasma flux tubes, they provide a semi-quantitative description of the 
ridge measurements (64,67,68,127). More detailed studies are feasible with real- 
istic hydrodynamic simulations (128). We note that the PHENIX collaboration 
which also observed ridge-like structures (129) have presented preliminary data 
showing that ridge-like structures disappear as the p± of the trigger particle is 
increased (130); this result is in qualitative agreement with expectations in the 
Glasma flux tube picture. For a recent critical evaluation of this and alternative 
models, we refer the reader to (131). 
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At the LHC, very low values of x will be probed in p+p, A+A (including the 
photo-production induced dynamics of peripheral A+A collisions) and possibly 
p+A collisions through a wide range of final states and diverse kinematic ranges. 
In nuclei, Q 2 will be large; estimates range from ~ 2.6-4 GeV 2 in central collisions 
to ~ 10 GeV 2 for y = ±3 units. The picture of strongly correlated albeit weakly 
coupled dynamics of glue in the CGC EFT outlined here will be tested as never 
before. Detailed tests of BFKL dynamics and possibly even the CGC are feasible 
in p+p collisions in studies respectively of Mueller-Navelet (132) and forward jets. 
Diffractive final states, while always challenging to interpret, offer opportunities 
to understand deeply how saturated gluons generate rapidity gaps. A strong test 
of CGC dynamics will already be available in "Day 1" physics of bulk dynamics 
in A+A collisions. More subtle tests of the RG flow of multi-parton correlators 
are available in a variety of final states ranging from quarkonia to electromagnetic 
probes to long range rapidity correlations. In p+A collisions at the LHC, many 
of the patterns seen in forward single inclusive and di-hadron correlation spectra 
will be present already at central rapidities and will be much more dramatic at 
forward rapidities. Finally, the role of Glasma dynamics relative to that of the 
Quark-Gluon Plasma in A+A collisions at the LHC is still unclear. If RHIC has 
reached the perfect hydrodynamic limit of maximal flow, what happens at the 
LHC? 

Studies of the dynamics of gluon saturation and the CGC that are comple- 
mentary to the LHC can be performed at a future Electron- Ion Collider (EIC) 10 . 
Firstly, with a wide lever arm in Q 2 , the dynamics of saturation can be studied 
with precision in the regime where Q 2 ~ Q 2 3> A 2 CD ; this is difficult to achieve 
at a hadron collider. Secondly, clarifying what aspects of the dynamics probed 
are universal and novel calls for an electron probe. The factorization theorems 
we have developed here suggest that the density matrices W are universal-how 
can we confirm this and cleanly extract their rich dynamics? These issues are not 
merely academic and are strongly reflected in the structure of final states. For 
example, rapidity gaps are a large fraction of the cross-section at HERA but are 
a much smaller contribution in hadronic collisions at Fermilab, demonstrating a 
breakdown of factorization (133). The physics case for an EIC/LHeC has been 

10 Current proposals include the EIC in the eRHIC version at BNL and the eLIC version at 

Jlab, the LHeC proposal at CERN and an electron-ion collider at FAIR in GSI. 



32 Frangois Gelis et al. 

outlined in (134) and in a number of white papers and reports. 

In summary, the prospects for unambiguous discovery and exploration of a 
novel many body QCD regime of gluon saturation are very bright in the next 
decade. The theoretical status of studies of this regime within the framework of 
the CGC EFT are increasingly robust albeit may questions remain. It is hoped 
that experiments, as usual, will provide definitive answers. 
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